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ABSTRACT 

Hubble Space Telescope (HST ) observations of the nearby (3.22 pc), K2 V star e Eridani have been combined 
with ground-based astrometric and radial velocity data to determine the mass of its known companion. We model 
the astrometric and radial velocity measurements simultaneously to obtain the parallax, proper motion, perturbation 
period, perturbation inclination, and perturbation size. Because of the long period of the companion, e Eri b, we 
extend our astrometric coverage to a total of 14.94 years (including the three year span of the HST data) by including 
lower-precision ground-based astrometry from the Allegheny Multichannel Astrometric Photometer. Radial velocities 
now span 1980.8 - 2006.3. We obtain a perturbation period, P = 6.85 ± 0.03 yr, semi-major axis a = 1.88 ± 0.20 
mas, and inclination i = 30?1 ± 3?8. This inclination is consistent with a previously measured dust disk inclination, 
suggesting coplanarity. Assuming a primary mass M* = 0.83M©, we obtain a companion mass M = 1.55 ± 0.24 Mj up . 
Given the relatively young age of e Eri (^800 Myr), this accurate exoplanet mass and orbit can usefully inform future 
direct imaging attempts. We predict the next periastron at 2007.3 with a total separation, p = 0'.'3 at position angle, 
p. a. = -27°. Orbit orientation and geometry dictate that e Eri b will appear brightest in reflected light very nearly 
at periastron. Radial velocities spanning over 25 years indicate an acceleration consistent with a Jupiter-mass object 
with a period in excess of 50 years, possibly responsible for one feature of the dust morphology, the inner cavity. 
Subject headings: astrometry — interferometry — stars: individual (e Eri) — stars: radial velocities 
— stars: late-type — stars: distances — extrasolar planets: masses 



1. INTRODUCTION 

e Eridani ( = HD 22049 = HIP 16537 = HR 1084 = 
PLX 742), with a spectral type of K2V, is one of the 
nearest solar-type stars with a distance of about 3.2 pc. 
It is slightly metal-poor (Fe/H = -0.13± 0.04, Santos 
et al. 2004, Laws et al. 2003). Its proximity makes it 
a prime target for future extrasolar planet direct imag- 
ing efforts. The success of these efforts will depend on 
knowing exactly where to look, requiring accurate orbital 
elements for the companion. It will depend on the mass 
of the planetary companion, and will depend on the age 
of the system. Younger and more massive gas giant plan- 
ets are predicted to be brighter (Hubbard et al. 2002). If 
young enough, the intrinsic luminosity of e Eri b might 
be greater than its brightness in reflected light. However, 
planetary mass objects with the age of e Eri change in- 
trinsic luminosity by a factor of nearly 100 between 1 
Mj up and 7 Mj up . Hence the need for a more precise 
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companion mass. 

e Eri has been the subject of multiple radial velocity 
(RV) planet searches. Walker et al. (1995), using mea- 
surements spanning 11 years, found evidence for « 10 yr 
variation with an amplitude of 15 m s" 1 . These results 
were substantiated by Nelson & Angel (1998) using an 
analysis of the same data set. Cumming et al. (1999) 
analyzed 11 years of RV data on this star taken at Lick 
Observatory and found significant variations with com- 
parable amplitude but with a shorter period of 6.9 years. 
Because of the high level of magnetic activity for e Eri 
(inferred from chromospheric activity), these RV vari- 
ations were largely interpreted as arising from a stel- 
lar activity cycle. The McDonald Observatory Planet 
Search Program (Cochran & Hatzes 1999) has monitored 
e Eri since late 1988. The McDonald results in combina- 
tion with these other surveys, along with data from ESO 
(Endl et al. 2002), confirmed the presence of long period 
RV variations and demonstrated that the most likely ex- 
planation for the observed RV variations was the pres- 
ence of a planetary companion with a period, P = 6.9 y. 
Details of this analysis are given in Hatzes et al. (2000). 

We obtained three years of astrometry with HST with 
millisecond of arc precision which we combined with ra- 
dial vel ocity data as we have in previous planeta r y mas s 
studies Benedict et al. I ()2002j) : iMcArthur et all (2004 ). 
Just as in the case of the Hatzes et al. (2000) radial 
velocity analysis, where less-precise data extended the 
observation span and allowed a companion detection, we 
anticipated that less-precise astrometry with a 14 year 
baseline from the Allegheny Multichannel Astrometric 
Photometer (MAP) astrometry project would improve 
the astrometric result. Gatewood (2000) reported in a 
meeting abstract an inclination, i = 46° ± 17° and a com- 
panion mass Mb — 1.2±0.3Mj up obtained with the MAP 



2 



data alone. In this study we combine the MAP data with 
the HST data, only to improve the determination of the 
proper motion of e Eri. The parallax and proper motion 
must be removed as accurately as possible to determine 
the perturbation orbit of e Eri, which, when combined 
with an estimate of the mass of e Eri, will provide the 
mass of the companion, e Eri b. 

This paper presents a mass of the planet orbiting 
e Eri discussed in Hatzes et al. (2000), not the far 
long er period object inferred fro m dust-disk morphol- 
ogy l)Quillen &: Thorndike 112002]) . Our mass is derived 
from combined astrometric and radial velocity data, con- 
tinuing a scries presenting accurate masses of planetary 
companions to n earby stars. Previous r esults include the 
m ass of Gl 876b (jBengdi ct et al. Il2002|) and of p 1 Cancri 
d (JMcArthur et al. 112004^) . 

In Section 2 we briefly review the astrometers and dis- 
cuss the data sets coming from each, and identify our 
many sources for radial velocities. In Section 3 we present 
the results of extensive spectrophotometry of the astro- 
metric reference stars, information required to correct 
relative parallax to absolute. In Section 4 we briefly dis- 
cuss our astrometric modeling and the quality of our re- 
sults as determined by residuals. In Section 5 we review 
our radial velocity data. In Section 6 we derive an abso- 
lute parallax and relative proper motion for e Eri, those 
nuisance parameters that must be removed to determine 
the perturbation orbital parameters. We finally estab- 
lish the perturbation orbital parameters and, combined 
with an estimate of the mass of e Eri, estimate a mass for 
e Eri b. We discuss system age, dust, and companion de- 
tectability in Section 7, and summarize our conclusions 
in Section 8. 

2. THE ASTROMETERS AND OBSERVATIONAL DATA 

2.1. HST FGSlr 

We used HST Fine Guidance Sensor lr (FGSlr) 
to carry out our space-based astrometric observations. 
Nelan et al. (2003) provides a detailed overview of FGSlr 
as a science instrument. Benedict et al (2002b) describe 
the FGS3 instrument's astrometric capabilities along 
with the data acquisition and reduction strategies used in 
the present study. We use FGSlr for the present study 
because it provides superior fringes fr om which to ob- 
tain t arget and reference star positions l|Mc Arthur et all 
2002). 

Table ^ presents a log of HST FGS observations. 
Epochs 2-4 contain multiple data sets acquired contigu- 
ously, the time span less than a day. Each time is that 
of the first observation within each epoch. Each distinct 
observation set typically contains five measurements of 
e Eri. The field was observed at multiple spacecraft roll 
values, and e Eri had to be placed in different, non-central 
locations within the FGSlr FOV to accommodate the 
less than optimal distribution of reference stars. FGS 
photometric sensitivity depends on location within the 
FOV {e.g. Benedict et al. 1998) and depends on time, 
as the FGS PMT age. Given the faintness of our ref- 
erence stars, we could not use them to provide a high- 
precision flat field. Hence, we could not extract milli- 
magnitude photometry (c.f. Benedict et al. 2000a) to 
monitor e Eri stellar activity. 

2.2. The Allegheny MAP 



The Multichannel Astrometric Photometer (MAP) 
and associated observation and reduction procedures are 
described by Gatewood (1987). The observational pro- 
gram utilizing the MAP began in 1986, but reluctance to 
observe at the low e Eri declination over the city of Pitts- 
burgh delayed initiation of its observation until January 
of 1989. Despite the reduced precision and rate of suc- 
cessful observation, the field remained on the MAP pro- 
gram until the installation of new instrumentation early 
in 2004 (Gatewood 2004). 

MAP observations of the brightest stars utilize either 
a specially filtered 12th channel (e.g. Gatewood & Han 
2006) or a divide by 16 feature on channel 2, which re- 
duces the count sufficiently for the 16 bit counters. With 
an R magnitude of 3.0 e Eri could be placed on channel 
2 without counter flooding. The other 10 channels were 
assigned to the MAP reference stars noted below. An ob- 
servation consists of 4 x 11 minute sweeps of the ruling 
across the field with probe rotations and ruling rotations 
to reduce systematic error (Gatewood 1987). Thus each 
observation consists of approximately 22 minutes of in- 
tegration on each axis. Table |21 presents a log of MAP 
observations. 

2.3. Radial Velocities 

The radial velocity data used include all data described 
in Hatzes et al. (2000), in addition to more recent data 
from McDonald Observatory and ESO. These data now 
span over 25 years. All sources are listed in Table |3J 
Briefly, McDonald Phases I, II, and III are all data ob- 
tained with the 2.7m Smith telescope. The phases cor- 
respond (I) to early velocities referenced to atmospheric 
O2, (II) velocities obtained with an I2 cell, and (III) ve- 
locities obtained with the McDonald 2d-Coude spectro- 
graph (Tull et al. 1995) and an I2 cell. The only new ra- 
dial velocity data included in this new study are from the 
McDonald Observatory Phase III 2.7m program. They 
are listed in Tabled Note that the errors associated with 
these data are larger than typicall y produced by this tele- 
scope/spectrograph combination ijEndl et al. 112006^1 due 
to high levels of e Eri stellar activity. 

3. e Eri ASTROMETRIC REFERENCE FRAMES 

Any prior knowledge concerning the fifteen stars in- 
cluded in our reference frame eventually enters our mod- 
eling as observations with error, and yields the most ac- 
curate parallax and proper motion for the prime target, 
e Eri. These periodic and non- periodic motions must be 
removed as accurately and precisely as possible to obtain 
the perturbation inclination and size caused by e Eri b. 

3.1. The MAP Reference Frame 

Figure ^ shows the distribution of the ten reference 
stars in the 130 MAP e Eri observation sets (Table |2J). 
Note that the areal coverage is approximately 0?6x0?6, 
allowing for the use of relatively bright reference stars, 
well-distributed around the prime target e Eri, in con- 
trast to the case for the HST FGS (below). 

3.2. The FGS Reference Frame 

Figure shows the distribution in FGSlr pickle coor- 
dinates of the 52 sets of five reference star measurements 
for the e Eri field. The arcing pattern is enforced by the 
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requirement that HST must roll to keep its solar panels 
fully illuminated throughout the year. To ensure access 
to all reference stars for every observation set, it was not 
possible to keep e Eri (symbol, x) located in the cen- 
ter of the FGSlr FOV. At each epoch we measured each 
reference stars 2-4 times, and e Eri 4-5 times. 

3.3. Absolute Parallaxes for the Reference Stars 

Because the parallax determined for e Eri is mea- 
sured with respect to reference frame stars which have 
their own parallaxes, we must either apply a statistically- 
derived correction from relative to absolute parallax 
l)van Altena Lee & Hoffleit"l l995. Yale Parallax Catalog, 
YPC95), or estimate the absolute parallaxes of the refer- 
ence frame stars. In principle, the colors, spectral type, 
and luminosity class of a star can be used to estimate the 
absolute magnitude, My, and V-band absorption, Ay. 
The absolute parallax is then simply, 



7Tabs = 10 



-{V-M v +b-A v )/b 



(1) 



3.3.1. Reference Star Photometry 



Our bandpasses for reference star photometry include: 
V (from FGSlr), and JHK from 2MASS 10 . The JHK 
values have been transformed to the Bessell & Brett 
(1988) system using the transformations provided in Car- 
penter (2001). Table lists VJHK photometry for the 
target and reference stars indicated in Figures H and 12 
Figure 01 contains a (J — K) vs. (V — K) color-color di- 
agram with reference stars and e Eri labeled. Schlegel 
et al. (1998) find an upper limit Ay^ 0.1 towards e Eri. 
In the following we adopt (Ay) — 0.0, but increase the 
error on reference star distance moduli by 0.1 magnitudes 
to account for absorption uncertainty. 

The derived absolute magnitudes are critically depen- 
dent on the assumed stellar luminosity, a parameter im- 
possible to obtain for all but the latest type stars using 
only Figure 13 To confirm the luminosity classes we ob- 
tain UCAC2 proper motions ijZacharias et al. Il2004|) for 
a one-degree-square field centered on e Eri, and then it- 
eratively employ the technique of reduced proper motion 
lYong fc Lambert Ejffll ICouTri fc Morgan 1200^ to dis- 
criminate between giants and dwarfs. The end result of 
this process is contained in Figure 

3.3.2. Adopted Reference Frame Absolute Parallaxes 

We derive absolute parallaxes using our estimated 
spectral types and luminosity class and My values from 
Cox (2000). Our adopted input errors for distance mod- 
uli, (m — M)o, are 0.4 mag for all reference stars (ex- 
cept ref-2 and -6, as discussed below). Contributions 
to the error are a small but undetermined Ay and er- 
rors in My due to uncertainties in color to spectral type 
mapping. We estimate a spectral type for reference star 
ref-6 only through its apparent magnitude, hence, the 
larger error in its distance modulus. Ref-2, which in 
Figure straddles the gap between giants and dwarfs, 
was finally typed K4V, because the x 2 (from modeling 
of the reference frame) significantly decreased, using that 
typing. Its input parallax error was also increased. All 

10 The Two Micron All Sky Survey is a joint project of the Uni- 
versity of Massachusetts and the Infrared Processing and Analysis 
Center / California Institute of Technology 



reference star absolute parallax estimates are listed in 
Table Individually, no reference star absolute par- 
allax is better determined than — = 18%. The aver- 
age input absolute parallax for the reference frame is 
(^abs) = 4.9 mas, a quantity known to ~ 5% (stan- 
dard deviation of the mean of fifteen reference stars). 
We compare this to the correction to absolute parallax 
discussed and presented in YPC95 (Sec. 3.2, Fig. 2). 
Entering YPC95, Fig. 2, with the Galactic latitude of 
e Eri, b = —48°, and average magnitude for the reference 
frame, (Kef) = 12.2, we obtain a correction to absolute 
of 2.3 mas, considerably different. Rather than apply a 
model-dependent correction to absolute parallax, we in- 
troduce our spectrophotomctrically-estimated reference 
star parallaxes into our reduction model as observations 
with error. 

4. THE ASTROMETPJC MODEL 

The e Eri reference frame contains fifteen stars and 
has been measured by two different astrometers, FGSlr 
and MAP. The only object in common is e Eri. From 
these positional measurements we determine the scale, 
rotation, and offset "plate constants" relative to an ar- 
bitrarily adopted constraint epoch for each observation 
set. As for all our previous astrometric analyses, we em- 
ploy GaussFit (Jefferys et al. 1988) to minimize \ 2 - The 
solved equations of condition for the e Eri field are: 

x =x + lcx(B - V) - AXFx (2) 

y' = y + lcy(B -V)- AXFy (3) 

£ = Ax' + By' + C - fi a At - P a n (4) 

r) = Dx' + Ey' + F - [i s At - P s tt (5) 

for FGSlr data and 

£ = Ax + By + C - P a n - fi a At (6) 

r) = Dx + Ey + F - P s ir - fx s At (7) 

for the MAP data. Identifying terms, x and y are the 
measured coordinates from HST and the MAP; (B — 
V) represents the (B — V) color of each star, estimated 
from its spectral type, Ay, and (J — K) color listed in 
Table IH1 and lex and Icy are the lateral color corrections, 
applied only to FGSlr data. Here AXFx and AXFy are 
the cross filter corrections in x and y, applied to the 
observations of e Eri in FGSlr. A, B, D and E are scale- 
and rotation plate constants, C and F are offsets; \i a and 
Us are proper motions; At is the epoch difference from 
the mean epoch; P a and Pg are parallax factors; and tt 
is the parallax. We obtain the parallax factors from a 
JPL Earth orbit predictor (Standish 1990), upgraded to 
version DE405. Orientation to the sky for the FGSlr 
data is obtained from ground-based astrometry (2MASS 
Catalog) with uncertainties of 0?01. 

4.1. Assessing Reference Frame Residuals 

Histograms of the MAP residuals (Figure |3J indicate 
per-observation precision of ~7 mas. Because we are 
seeking the signature of a perturbation over three times 
smaller than that per-observation precision, the MAP 
data were only used to lower the errors on parallax and 
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proper motion, not to establish any perturbation param- 
eters. As for the F GS data, the Optical Fie ld Angle Dis- 
tortion calibration McArth ur et al. I (12002^ reduces as- 
built HST telescope and FGSlr distortions with mag- 
nitude ~ 1" to below 2 mas over much of the FGSlr 
field of regard. From histograms of the FGS astrometric 
residuals (Figure EJ) we conclude that we have obtained 
correction at the ~ 1 mas level. The reference frame 
'catalogs' for MAP and FGSlr in £ and r\ standard co- 
ordinates (Table IHJ) were determined with < cr^ >= 1.0 
and < a v >= 1.3 mas (MAP), and < er^ >= 0.3 and 
< cr v >= 0.2 mas (FGS). 

5. RADIAL VELOCITIES 

Measurements from four planet search groups were in- 
cluded in our modeling. Table lists the source, cov- 
erage, technique, number of observations and the rms 
deviation from the final orbit of these observations. The 
weighting of the RV data was carefully evaluated with in- 
dependent modeling and significant outliers were filtered. 
For example - if five data points were taken in succession, 
all assigned with the same weight and one point was 50 
m s _1 offset from the others, that point was discarded 
for this solution. Initially these points were merely re- 
weighted, but later discarded as spurious data. The total 
number of observations so discarded was less than 1% of 
the aggregate data. This improved our goodness of fit 
(x 2 /degrees of freedom) measurement of the modelling 
for the RV data set from 0.94 in the announcement paper 
to 0.30 in the current analysis. 

6. e Eri PARALLAX, PROPER MOTION, AND 
PERTURBATION ORBIT FROM ASTROMETRY AND 
RADIAL VELOCITIES 

Solving for relative parallax, proper motion, and or- 
bital motion, using astrometry and radial velocities si- 
multaneously, the model now becomes, 

£ = aX + bY + c-P x *w- fj, x *t- ORBIT x (8) 
77 = -bX + aY + / - P y *7r - Hy *t - ORBITy (9) 

where ORBIT is a function (through Thiele-Innes con- 
stants) of the traditional astrometric and radial velocity 
orbital elements listed in Table UTI 

The period (P), the epoch of passage through perias- 
tron in years (T), the eccentricity (e) and the angle in the 
plane of the true orbit between the line of nodes and the 
major axis (u>), are constrained to be equal for the ra- 
dial velocity and astrometry portions of the model. Only 
radial velocity provides information with which to deter- 
mine the half- amplitudes (K\) and 7, the systemic ve- 
locity. Combining radial velocity observations from dif- 
ferent sources is possible with GaussFit, which has the 
ability to simultaneously solve for many separate velocity 
offsets (because velocities from different sources are rel- 
ative, having differing zero points), along with the other 
orbital parameters. 

We force a relationship between t he astrometry and the 
radial velocity by a constraint from lPourbaix fc Jorissenl 
(2000) 

a A sin 1 = PK^l-e 2 ) 1 ' 2 

ir abs 2tt x 4.7405 1 ' 



where quantities derived only from astrometry (parallax, 
TTabs, primary perturbation orbit size, a a, and inclina- 
tion, i) are on the left, and quantities derivable from 
both (the period, P and eccentricity, e), or radial veloc- 
ities only (the radial velocity amplitude for the primary, 
Ki), are on the right. 

Combining RV measurements complete through 2006.3 
(Table |3J), all the astrometric measurements, and the 
Equation 10 constraint, we solve for parallax, proper 
motion, and the semi-major axis, orbit orientation, and 
orbit inclination for the perturbation caused by the com- 
panion. For the parameters critical in determining the 
mass of e Eri we find a parallax, 7r a (, s = 311.37 ± 0.10 
mas and a proper motion 976.54 ± 0.1 mas y _1 in posi- 
tion angle 269?0 ± 0?6. Table ITUI compares values for the 
parallax and proper motion of e Eri from HST and HIP- 
PARCOS. We note satisfactory agreement. Our precision 
and extended study duration have significantly improved 
the accuracy and precision of the parallax and proper 
motion of e Eri. 

At this stage we can assess the reality of any e Eri per- 
turbation by plotting residuals to a model that does not 
include an orbit. Figure shows the X and Y compo- 
nents of only the higher-precision astrometry FGS resid- 
uals plotted as grey dots. The lower precision MAP data 
were not considered in the determination of the orbital 
parameters. We also plot normal points formed from 
those dots at nine epochs. Finally, each plot contains as 
a dashed line the X and Y components of the perturba- 
tion we find by including an orbit in our modeling. 

We find a perturbation size, cxa = 1-88 ± 0.19 mas, 
and an inclination, i= 30?1 ± 3?8. These, and the other 
orbital elements for the perturbation, are listed in Ta- 
ble El with 1-er errors. Errors generated by GaussFit 
(Jefferys et al. 1988) come from a maximum likelihood 
estimation that is an approximation to a Bayesian max- 
imum a posteriori estimator with a flat prior (Jefferys 
1990). Figure [S] illustrates the Pourbaix and Jorrisen re- 
lation (Equation 10) between parameters obtained from 
astrometry (left-side) and radial velocities (right side) 
and our final estimates for a a and i. As seen in Ta- 
bles ED and 111! most of the errors of the terms in Equa- 
tion 10 are quite small. In essence, our simultaneous so- 
lution uses the Figure |H1 curve as a quasi-Bayesian prior, 
sliding along it until the astrometric and radial veloc- 
ity residuals are minimized. Gross deviations from the 
curve are minimized by the high precision of many of 
the terms in Equation 10. Figure El contains all radial 
velocity measures and the predicted velocity curve from 
the simultaneous solution. Compared to the typical per- 
turbation radial velocity curve (e.g. Hatzes et al. 2005, 
McArthur et al. 2004, Cochran et al. 2004), Figure|H|ex- 
hibits far more scatter about the derived orbit. There 
are two reasons for this. The perturbation amplitude is 
small (Kl=18.5 m s _1 ), and e Eri is an active star, as 
discussed in Hatzes et al. (2000). Reiterating their con- 
clusions, none of the activity cycles have periods com- 
mensurate with the planetary perturbation period. Fig- 
ure n21 presents the astrometric residuals and the derived 
perturbation orbit for the primary star, e Eri. Stellar 
activity has even less of an effect on astrometry at our 
level of precision. A star spot covering 30% of the sur- 
face would induce a photocenter shift of less than 0.2 mas 
(Sozzetti 2005). The astrometry confirms the existence 
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of the companion. 

Our analysis of the radial velocities (now spanning over 
25 years, all shown in Figure |3J) included a linear drift 
term, a change in velocity as a function of time. This 
drift is clearly seen in the overplotted final radial ve- 
locity curve, and amounts to 0.32 ± 0.05 m s^ 1 yr" 1 . 
Such a change can be caused by longer-period cqmpan - 
ions and/or secular acceleration ijKurster k, Endl 1 2004). 
The secular acceleration expected for e Eri is (van de 
Kamp, 1967) quite small, 0.07 m s _1 yr _1 . The trend 
we find is over four times larger than the predicted sec- 
ular acceleration, and is approximately the acceleration 
one might find for a planet similar in mass to e Eri b, 
but with a 50-100 year period. Figure 9 shows sufficient 
enough overlap among the many velocity data sets that 
the trend is unlikely an artifact due to a mismatch in the 
center of mass velocity offsets (discussed above) obtained 
for each set. Typical offset random error is ~ 1 m s _1 . 
While this acceleration is not a detection of the longer- 
period companion (40 < a < 60 AU) invoked by Quillan 
& Thorndike (2002) and Ozernoy et al. (2000) to mod- 
ify the dust distribution as discussed below in Section 
7, it may (with a semimajor axis 10-20 AU) be at least 
partially responsible for the inner cavity in the dust disk 
distribution imaged by Greaves et al. (2005). The astro- 
metric motion over 15 years due to this possible tertiary 
would be of order 3 mas and difficul t to se parate from 
proper motion, e.g., IBlack fc Scareiel flQSl . 

The planetary mass depends on the mass of the pri- 
mary sta r, for which we have adopted M* = 0.83 ± 
0.05M Q Ipi Folco et aJ. 1 12004^ . For this A/* we find 
Mf, = 1.55 ± 0.24M,/ up . The companion is clearly an 
extrasolar giant planet. In Table IT2l the mass value, 
incorporates the present uncertainty in M* . Until e Eri b 
is directly detected, its radius is unkn own. Fro m a review 
of exoplanet masses and radii (Guillot 2005), a radius 
of R — lRjup seems reasonable. 

Our eccentricity value, e = 0.70, allows for a significant 
difference in separation between star and exoplanet at 
apastron compared to periastron. At time of periastron 
passage, T = 2007.29, we predict a separation ft'3 ± 
O'.'l at a position angle of -27°. At the next apastron, 
to occur 2010.71, the separation should be l'/8 ± 0'.'4 at 
position angle of 153°. The dominant sources of error for 
the separations are the eccentricity (6%) and the e Eri b 
planet mass (15%). 

7. DISCUSSION 

Our accurate mass and orbital parameters for this 
planetary companion have value for future direct imag- 
ing projects. We now know where near e Eri to look for 
e Eri b. We would now like to know when to look, what 
bandpass is best, and what we can expect to see. As 
stated previously system age, companion mass, and or- 
bital geometry are critical parameters when estimating 
visibility. 

A high level of chromospheric activity is seen for 
e Eri (e.g. Gray & Baliunas 1995), and is consistent with 
a relatively young age; < 1 Gyr (Soderblom & Dappen 
1989). Saffe et al. (2005) used the calibrations of Don- 
ahue (1993) and Rocha-Pinto & Maciel (1998) (which 
corrected the age with an effect from stellar metallic- 
ity) to estimate ages of 0.66 and 0.82 Gyr, respectively. 
Henry (1986) derived from Call lines a value of 0.8 Gyr. 



Song et al. (2000) used Li abundances with its position 
in the H-R diagram and kinematics to derive a value of 
0.73 ± 0.2 Gyr. Di Folco et al. (2004) estimated the 
age at 0.85 Gyr, a value obtained through the measure- 
ment of the radius of e Eri by long-baseline interferome- 
try. Their modeling is consistent with a primary mass of 
M* = 0.83 ± 0.05Af Q , an estimate that weakly depends 
on measured metallicity, which ranges -0.13 < [Fe/H] < 
-0.06 in the literature. 

Hubbard et al. (2002) predict the intrinsic luminosity 
of extrasolar giant planets as a function of mass and age. 
From their Figure 11, an age of 800 Mv lDi Folco et all 
(2004), and our planetary mass, M b = 1.55 ± 0.24M Jup , 
we find for e Eri b, L/L Q = 1.6 X 10 -8 . Using the Di Folco 
et al. (2004) T e// =5135 K, their radius, i?» = 0.743i? Q , 
our parallax, n a f, s = 311.37 mas, and a bolometric cor- 
rection, B.C. = -0.27, from Flower (1996), we find a dif- 
ference in bolometric magnitude of e Eri compared to 
the Sun of AMb i = +1.17. Hence, neglecting reflected 
light and orbital phase, e Eri b is 4.67xl0 -8 fainter in 
bolometric luminosity than e Eri. 

Sudarsky et al. (2005), Dyudina et al. (2005), and Bur- 
rows et al. (2004) discuss exoplanet apparent brightness 
in reflected host star light as functions of orbit geome- 
try, orbital phase, and cloud cover. Burrows et al. (2004) 
predict the full spectrum of e Eri b from 0.5 to 6/im, as- 
serting that the planet is too young for its atmosphere 
to contain condensed ammonia clouds. However, e Eri b 
should exhibit H2O clouds. They predict a maximum 
planet/host star flux ratio, \og w (F p i anet / F star ) ~ -7, at 
~ 4.5/im with a secondary peak, \og w (F p i anet / F star ) ~ 
—8, at ~ 1/im. Dyudina et al. (2005) predict that for 
oj = 30°, inclination, i = 30°, e = 0.5, and a Jupiter- 
like atmosphere, the planet /host star flux ratio is largest 
very shortly after periastron, late 2007. However, the 
separation remains small (^0'.'3). The inclination of the 
e Eri system, i = 30°, is likely to decrease the flux ratio 
by approximately a factor of two (Sudarsky et al. 2005), 
compared to a i = 90° edge-on orientation. Given the 
orientation of the orbit of e Eri b (its ascending node, W 
= 254°), the disk of e Eri b is most fully illuminated at 
apastron, but is three times further away from its pri- 
mary. 

The dusty rings or debris disks surrounding e Eri also 
suggest relative youth for the system. Photometric 
measurements from the IRAS satellite (Aumann, 1988) 
provided the first hint of dust around e Eri. Subse- 
quently, Submillimcter Common-User Bolometric Array 
(SCUBA) measurements were made between 1997 and 
2002. These measurements determined that the dust, 
distributed in a ring , is located 65 AU from the star 
lIGreaves et al. lll9M 12005). The sub -mm bolometer, 
SIMBA, provided observational confirmation of this ex- 
tended dust disk (Schutz et al. 2004). The STIS CCD 
camera on -ffSTtook deep optical images around e Eri in 
an effort to find an optical counterpart for the sub- 
millimeter observations. These measurements did not 
provide clear evidence for the detection of that optical 
counterpart, but did place a limit on the optical surface 
brightness of the dust, that it could not be brighter than 
approximately 25 STMAG arcsec -2 , which places con- 
straints on the nature and amount of the smallest dust 
grains (Proffit et al. 2004). 
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Observational and theoretical searches for the signa- 
ture of planetary /brown dwarf objects in the structure 
of the dust disk around e Eri are underway. Clumps 
seen in the ring are thought to come from the inter- 
action between the disk and a massive planetary body 
(Holland et al. 2003) . Adaptive optics on the Keck Tele- 
scope were used to search for extrasolar planets. These 
studies found no evidence of brown dwarf or planetary 
companions down to 5 Jupiter masses at the angular sep- 
arations comparable to that of the dust rings (Macin- 
tosh et al. 2003). Spitzer Space Telescope (SST) obser- 
vations made with the Multi-band Imaging Photometer 
(MIPS) and the InfrRed Spectograph (IRS) have con- 
firmed the disk and provided evidence for asymmetries 
in the structure of the disk that may have been caused by 
the gravitational perturbation of sub-stellar companions 
(Marengo et al. 2004). 

Two recent studies suggested that debris disks and 
long-period planets co-exist, with planetary bodies 
'sculpting' the disk, e Eri is the prototypical system. 
First, high-resolution modeling of the structure of the 
disk around e Eri predict an angular motion of the 
asymmetry of the disk of about 0?6-0?8 yr _1 (Ozervnoy 
et al. 2000). Secondly, Quillen & Thorndike (2002) car- 
ried out numerical simulations of dust particles captured 
in mean motions resonances with a hypothetical planet 
(e = 0.3 , M = 10~ 4 M , a = 40 AU) at periastron. 
These produced a dust distribution that agreed with the 
morphology of the dust ring around e Eri presented by 
Greaves et al. (1998, 2005). An investigation into the 
dynamics of the dust ring around e Eri (Moran, Kuchncr 
& Holman, 2004) concluded that the eccentricity of the 
dust released in the inner ring (<20 AU) could reveal 
patterns in the dust which could confirm the existence of 
the planet reported by Hatzes et al. (2000). 

We determined an inclination of i=30°l ± 3?2 for 
e Eri b. Our measured inclination is consistent with 
the previously measu red dust disk in c lination from 450 
and 850 /im maps of iGreaves et al. I l|1998l l2005|) . i = 
25°. This suggests that the dust disk and plane of 
the orbit of e Eri b are coincident and that the dust 
distribution is nearly circular. This provides support 
for hiearchica l accre tion models for planet formation 
iPollack et al. I l|1996ft . where coplanar dust and a debris 
disk are expected remnants of planet formation (Tsiganis 
et al. 2005). Lastly, e Eri b and the possible tertiary de- 
duced from the linear trend in the radial velocities (Sec- 
tion 6) would most likely eject particles that would spiral 
inward, and recent SCUBA submillimeter observations 
have shown that the center of the disk is relatively exca- 
vated of dust, with half or less of the signals seen in the 
ring (Greaves et al. 2005). 

8. CONCLUSIONS 

Analyzing three years of HST FGS and over 14 years 
of Allegheny Observatory MAP astrometry, we find an 
independently determined parallax and proper motion 
for e Eri that agree within the errors with HIPPARCOS. 

Astrometric observations with HST FGS, combined 
with long-duration Allegheny MAP astrometry and 
ground based radial velocities, have confirmed the ex- 
istence of the planet orbiting e Eri, first suggested by 
Walker et al. (1995), noted by Cummings et al. (1999), 
and finally announced by Hatzes et al. (2000). 



Combining the astrometry with radial velocities from 
six different sources, spanning 25 years, and applying 
the Pourbaix & Jorrisen constraint between astrome- 
try and radial velocities, we obtain for the perturb- 
ing object e Eri b a period, P=6.85 ± 0.02 y, inclina- 
tion, i=30?l ± 3?2, and perturbation semimajor axis, 
ola = 1-88 ± 0.19 mas. Assuming for e Eri a stellar mass 
M* = 0.83±0.05Af©, we obtain a mass for e Eri b, M& 
= 1.55±0.24Af/ up . This companion inclination matches 
the disk inclination determined by Greaves et al. (2005). 

Our astrometry predicts for e Eri b periastron passage 
at To = 2007.29, with a separation ~0'/3 in position an- 
gle -27°(a = 3.39 AU). At the next apastron, to occur 
2010.71, the separation should be 1'.'7 in position angle 
153°. The orbital geometry suggests that 2007.97 (late 
December 2007) is the most favorable time for direct de- 
tection in reflected light. For an e Eri age ~850 My 
and our determined mass, e Eri b will have an intrinsic 
luminosity L/L e = 1.6 x 10~ 8 , 4.67xl0~ 8 times fainter 
in bolometric luminosity than e Eri. 

Radial velocities spanning 25 years indicate a long- 
term linear trend, an acceleration consistent with a 
Jupiter-mass object with a period of 50-100 years. This 
is a possible detection of a tertiary companion respon- 
sible for a major feature of the dust morphology, the 
central cavity. 
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TABLE 1 
Log of e Eri FGS Observations 
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a MJD = JD - 2400000.5 

b Spacecraft roll as defined in Chapter 2, FGS Instrument Handbook Nclan ct aJ. (2003) 
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Log of e Eri MAP Observations 
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56 
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121 


52925.3774 


2003.7806 


57 


52227.2226 


2001.8692 


122 


52937.2087 


2003.8130 


58 


52265.0691 


2001.9728 


123 


52937.2413 


2003.8131 


59 


52893.3497 


2003.6929 


124 


52946.224 


2003.8377 


60 


52924.2726 


2003.7776 


125 


52946.2566 


2003.8378 


61 


52925.2747 


2003.7804 


126 


52946.2899 


2003.8379 


62 


52925.3420 


2003.7805 


127 


53000.058 


2003.9851 


63 


52937.2087 


2003.8130 


128 


53001.076 


2003.9879 


64 


52946.2240 


2003.8377 


129 


53001.1177 


2003.9880 


65 


53000.0580 


2003.9851 


130 


53001.1448 


2003.9881 



a MJD = JD - 2400000.5 
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TABLE 3 
The Radial Velocity Data Sets 



Data Set 


Coverage 


Technique 


N 


RMS 




(yr) 






(m s- 1 ) 


CFHT 


1980.81-1991.88 


HF cell 


48 


10.5 


Lick 


1987.69-1998.99 


Iodine cell 


54 


11.5 


McD 01 


1988.74-1994.81 


Telluric 


27 


15.2 


McD 4>Il 


1990.78-1998.07 


Iodine cell 


42 


11.7 


ESO 


1992.84-1998.02 


Iodine cell 


36 


9.6 


McD 0III 


1998.69-2004.86 


Iodine cell 


28 


7.4 






total 


235 





TABLE 4 



New McDonald </>III Radial Velocities 


mJD 




mJD 






(m s- 1 ) 




(m s- 1 ) 


51066.4339 


-3.6±6.4 


52539.4226 


-16.4±6.4 


51212.1671 


13.6 5.0 


52576.4378 


2.0 5.5 


51239.1133 


1.7 8.6 


52661.1445 


-11.8 5.2 


51449.4333 


6.5 5.0 


52931.3912 


-13.9 6.4 


51503.3574 


14.7 5.1 


52958.2449 


-0.4 5.4 


51529.1986 


-3.5 7.4 


52958.2481 


0.0 6.0 


51555.1645 


-0.2 5.2 


53016.2424 


2.9 6.0 


51775.4643 


1.3 5.9 


53016.2456 


2.1 6.0 


51809.4033 


-2.7 4.8 


53035.1757 


-5.2 5.6 


51917.2260 


2.9 5.5 


53075.0940 


3.0 5.7 


51984.0740 


-2.8 5.0 


53318.3124 


-10.1 6.6 


52142.4227 


-0.7 6.2 


53632.4520 


-1.2 5.3 


52142.4264 


-0.4 6.1 


53632.4550 


0.4 5.5 


52248.2921 


-14.0 5.9 


53689.3900 


-3.0 5.2 


52303.1266 


-16.1 7.7 


53745.2110 


6.3 5.1 


52328.1208 


2.1 5.6 


53809.0850 


9.0 5.5 


52330.1064 


-0.9 5.7 
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TABLE 5 

FGSlR AND MAP ASTROMETRIC REFERENCE STARS 



ID 


Catalog 


RA a (2000.0) 


Dec a 


V 


2MASS 


1 


e Eri 


53.232961 


-9 458295 


3.82 


03325591 -0927298 


2 




53 325845 


-9 467569 


15.61 


03331 820-0928032 


3 




53 31 2638 


-9.434314 


15.58 


03331 503-0926035 


4 




53.295306 


-9.421752 


16.12 


03331087-0925183 


5 




53 269681 


-9.418811 


16.41 


03330472-0925077 


6 




53.207917 


-9.449667 


16.48 




7 


GEN# +6.10280864 


53.374062 


-9.504524 


11.36 


03332977-0930162 


8 


HD 22130 


53.401866 


-9.349848 


9.52 


03333644-0920594 


9 


BD-09 696 


53.164328 


-9.267778 


10.69 


03323943-0916040 


10 


GEN# +6.10280861 


53.032108 


-9.514245 


11.60 


03320770-0930512 


11 


BD-10 699 


53.399013 


-9.586081 


10.24 


03333576-0935098 


12 


HD 21951 


53.001248 


-9.385150 


9.69 


03320029-0923065 


13 


BD-10 695 


53.320908 


-9.696887 


9.78 


03331701-0941487 


14 


BD-10 700 


53.466004 


-9.641330 


9.94 


03335184-0938287 


15 


BD-09 699 


53.429559 


-9.232722 


11.20 


03334309-0913577 


16 


2MASS 03320556-0945292 


53.023188 


-9.758118 


11.15 


03320556-0945292 



a Positions from 2MASS, except ID 6 from applying STScI Visual Target Tuner to the Digital Sky Survey. 
b Magnitudcs from FGSlr (ID 2-6), SIMBAD (ID 1), or MAP (ID 7- 16). 



TABLE 6 
V and Near-IR Photometry 



ID 


V 


K 


(J-H) 


(j -JO 


(V-K) 


1 


3.82+0.01 


1.82+0.05 


0.40+0.28 


0.48+0.07 


2.00+0.05 


2 


15.61+0.03 


13.01+0.03 


0.58+0.04 


0.72+0.04 


2.60+0.04 


3 


15.58+0.03 


14.09+0.06 


0.41+0.05 


0.44+0.07 


1.50+0.06 


4 


16.12+0.03 


13.88+0.05 


0.51+0.04 


0.59+0.06 


2.24+0.06 


5 


16.41+0.03 


14.21+0.08 


0.46+0.07 


0.62+0.09 


2.20+0.08 


6 


16.48 










7 


11.36+0.03 


9.93+0.02 


0.29+0.03 


0.33+0.03 


1.43+0.04 


8 


9.52+0.01 


8.63+0.02 


0.21+0.06 


0.24+0.03 


0.89+0.02 


9 


10.69+0.03 


8.85+0.02 


0.48+0.04 


0.54+0.03 


1.84+0.04 


10 


11.60+0.03 


10.43+0.02 


0.33+0.03 


0.41+0.03 


1.17+0.0 4 


11 


10.24+0.03 


8.73+0.02 


0.40+0.06 


0.43+0.04 


1.51+0.04 


12 


9.69+0.03 


8.98+0.02 


0.13+0.03 


0.18+0.03 


0.71+0.04 


13 


9.78+0.01 


7.87+0.02 


0.47+0.03 


0.55+0.02 


1.91+0.02 


14 


9.94+0.03 


7.69+0.03 


0.57+0.06 


0.66+0.04 


2.25+0.04 


15 


11.20+0.03 


8.13+0.03 


0.62+0.05 


0.74+0.03 


3.07+0.04 


16 


11.15+0.03 


9.43+0.02 


0.45+0.03 


0.49+0.03 


1.72+0.04 



TABLE 7 

ASTROMETRIC REFERENCE STAR ADOPTED SPECTROPHOTOMETRIC PARALLAXES 



ID 


Sp. T. a 


V 


M v 


m-M 


Tabs (mas) 


2 


K4V 


15.6 


7.1 


8.5+1 


1.9+1.0 


3 


G8V 


15.6 


5.6 


10.0+0.4 


1.0+0.2 


4 


K2V 


16.1 


6.5 


9.6+0.4 


1.2+0.2 


5 


K2V 


16.4 


6.5 


9.9+0.4 


1.0+0.2 


6 


K2V 


16.4 


6.5 


9.9+2 


1.0+0.9 


7 


GOV 


11.4 


4.4 


7.0+0.4 


4.1+0.7 


8 


F5V 


9.5 


3.5 


6.0+0.4 


6.3+1.2 


9 


K0V 


10.7 


5.9 


4.8+0.4 


11.0+2.0 


10 


G5V 


11.6 


5.1 


6.5+0.4 


5.0+0.9 


11 


G8V 


10.2 


5.6 


4.6+0.4 


11.8+2.2 


12 


F0V 


9.7 


2.7 


7.0+0.4 


4.0+0.7 


13 


K0V 


9.8 


5.9 


3.9+0.4 


16.7+3.1 


14 


K0III 


9.9 


0.7 


9.2+0.4 


1.4+0.3 


15 


K2III 


11.2 


2.7 


8.5+0.4 


2.0+0.4 


16 


G8V 


11.2 


5.6 


5.6+0.4 


7.8+1.4 



a Spectral types and luminosity class estimated from colors and reduced proper motion diagram. 
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TABLE 8 

e Eri and Reference Star Relative Positions 



ID 




V 


e 




rf 


l b 




3.73 


297. 0383±0. 0001 




35.6893±0.0001 


2 




15.57 


-25.6479±0.0002 




126.1290±0.0001 


3 C 




15.57 


0.0000±0.0002 




0.0000±0.0002 


4 




16.09 


52.8710±0.0003 




-55.1822±0.0002 


5 




16.40 


433. 1267±0. 0003 




-20.0260±0.0003 


6 




16.37 


140.6876±0.0004 




-81.2800±0.0002 


l d 




3.73 


715. 2729±0. 0002 




201. 7543±0. 0002 


7 




11.36 


1214.0919±0.0013 




35.0491±0.0014 


8 




9.52 


1313.1141±0.0005 




591.8986±0.0007 


9 




10.69 


469.2508±0.0011 




887.5650±0.0014 


10° 




11.60 


0.0000±0.0021 




0.0000±0.0025 


11 




10.24 


1302.5205±0.0007 




-258.4640±0.0011 


12 




9.69 


-109.9029±0.0008 




464.7282±0.0008 


13 




9.78 






nco/i-i-n nnri7 


14 




9.94 


1540.2424±0.0005 




-457.4732±0.0006 


15 




11.2 


1411.6967±0.0011 




1013.5480±0.0018 


16 




11.15 


-31.1988±0.0015 




-877.8886±0.0023 


a £ and rj are relative positions in arcseconds 








b cpoch 2002.614, J2000 










C RA 


= 53.312638, Dec = 


-9.434314, J2000 








d cpoch 1996.761, J2000 










e RA 


= 53.032108, Dec = 


-9.514245, J2000 














TABLE 9 










Reference Star Proper Motions 










Input (UCAC2) 




Final (HST) 




ID 


V 








M<s a 


2 


15.57 


0.0053±0.0078 


-0.0112±0.008 


0.0110±0.0002 


- 0.0076±0.0002 


3 


15.57 


0.0166±0.0076 


-0.0014±0.0076 


0.0072±0.0003 


- 0.0026±0.0002 


4 


16.09 


0.0044±0.009 


-0.0028±0.0079 


0.0080±0.0004 


- 0.0022±0.0002 


5 


16.4 


-0.0004±0.004 


-0.0044±0.004 


0.0070±0.0004 


0.0082±0.0002 


6 


16.37 


0.0064±0.003 


0.0083±0.003 


0.0014±0.0004 


-0.0030±0. 0003 


7 


11.36 


0.0253±0.005 


0.0117±0.005 


0.0221±0.0002 


0.0083±0. 0003 


8 


9.52 


0.0104±0.005 


-0.0025±0.005 


0.0088±0.0001 


- 0.0052±0.0001 


9 


10.69 


0.003±0.005 


-0.0178±0.005 


0.0035±0.0002 


- 0.0184±0.0003 


10 


11.6 


-0.0015±0.005 


-0.0074±0.005 


-0.0020±0.0004 


- 0.0029±0.0004 


11 


10.24 


0.0101±0.005 


-0.0094±0.005 


0.0126±0.0001 


- 0.0089±0.0002 


12 


9.69 


0.0095±0.005 


-0.0038±0.005 


0.0109±0.0001 


- 0.0074±0.0002 


13 


9.78 


0.0549±0.005 


-0.037±0.005 


0.0579±0.0001 


- 0.0326±0.0001 


14 


9.94 


-0.0003±0.005 


-0.0079±0.005 


-0.0017±0.0001 


- 0.0102±0.0001 


15 


11.2 


0.0117±0.005 


-0.0063±0.005 


0.0108±0.0002 


- 0.0042±0.0003 


16 


11.15 


0.0073±0.005 


-0.0136±0.005 


0.0052±0.0003 


- 0.0158±0.0004 



a /i a and fig are relative motions in arcsec yr 



TABLE 10 
e Eri Parallax and Proper Motion 



Parameter 


HST 


MAP 


Combined 


Study duration 


2.92 y 


14.94 y 




number of observation sets 


46 


130 




reference star (V) 


16.0 


10.2 


12.2 


reference star ({B — V)) 






0.9 a 


Absolute Parallax 15 






311.37 ± 0.11 mas 


Relative Proper Motion 






976.54 ± 0.1 mas y" 1 


in pos. angle 






269? ± 0?6 


HIPPARCOS Absolute Parallax 
HIPPARCOS Proper Motion 
in pos. angle 






310.74 ± 0.85 mas 
976.52 ± 1.9 mas y- 1 
271? 1 ± 3? 8 



a Estimated from VJHK photometry with Ay = 0.0. 

b Value from modeling RV and HST and MAP astrometry simultaneously 



13 



TABLE 11 

Orbital Elements of e Eai Perturbation Due to e Eri b 



Parameter Value 

a A 1-88 ± 0.20 mas 

a A sini 3.02e-3 ± 0.32e-3 AU 

P 2502 ± 10 d 

P 6.85 ± 0.03 yr 

T 54207 ±7 mJD 

T 2007.29 ± 0.02 y 

e 0.702 ± 0.039 

i 30? 1 ± 3?8 

Q 74° ± 7° 

w 47° ± 3° 

Ki 18.5 ±0.2 m s" 1 

M* 0.83±0.05M© 



TABLE 12 

e Eri b Parameters 



Parameter Value 

a (AU) 3.39 ± 0.36 

ST 254° 

u>' 47° 
Mass Function (M ) 5.9e-10 ± 1.0e-10 

Msinj (Mj UP ) a 0.78 ± 0.08 

M (M JUP ) b 1.55 ± 0.22 

M (M JUP ) c 1.55 ±0.24 



a derived from radial velocity alone 

b derived from radial velocity and astromctry, using Msini/sini 
c derived from radial velocity and astrometry, using m\l(ra\ +m,2) 2 = a 3 /P 2 ; 
includes host star mass uncertainty. 



14 



> 



1000 - 




e Eri & MAP Reference Stars 






9 


15 


500 - 


12 




8 


- 




1 






10 




7 


-500 - 






11 

14 


1000 - 


16 




13 



-1000 -500 500 

X (arcsec) 



Fig. 1. — e Eri and MAP reference stars on the sky. Each star is identified by the number listed in Table 151 The inner box indicates the 
FGSlr reference frame coverage for epoch 2 in Tabic HI 
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Fig. 2. — HST e Eri and reference frame observations in FGSlr pickle coordinates. The symbol shape identifies each star listed in 
Tablc|S] Note that the position of e Eri (x) within the FGSlr FOV is not fixed at the center. 
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Fig. 3. — (J — K) vs. (V — K) color-color diagram for stars identified in Table 151 The dashed line is the locus of dwarf (luminosity class 
V) stars of various spectral types; the dot-dashed line is for giants (luminosity c lass III), T he reddening vector indicates Ay = 1.0 for the 
plotted color systems. Along this line of sight maximum extinction is Ay^ 0.1 ISchlegel et ai. 1119981) . 
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Fig. 4. — Reduced proper motion diagram for 9,041 stars in a 6° field centered on e Eri. Star identifications are in Table 151 For a given 
spectral type, giants and sub-giants have more negative Hk values and are redder than dwarfs in (J — K). Hk values are derived from 
'Final' proper motions in Table l9l The small cross at the lower left represents a typical (J — K) error of 0.04 mag and Hk error of 0.17 
mag. The horizontal dashed line is a giant-dwarf demarcation derived from a statistical analysis of the Tycho input catalog (Ciardi 2004, 
private communication). Ref-14 and -15 are likely luminosity class III. 
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Fig. 5. — Histograms of x and y residuals obtained from modeling the MAP observations of e Eri and the MAP reference frame with 
equations 6 and 7. Distributions are fit with Gaussian distributions. 
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Fig. 6. — Histograms of x and y residuals obtained from modeling the FGS observations of e Eri and the FGS reference frame with 
equations 4 and 5. Distributions are fit with Gaussian distributionss. 
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Fig. 7. — X and Y componen ts o f the perturbation orbit for e Eri as a function of time. The dashed line is the orbit described by the 
orbital elements found in Table 1111 The dots are all the e Eri observation residuals to a model that does not contain orbital motion. 
The circles are normal points formed from the individual observation residuals, where the circle size is proportional to the number of 
observations forming the normal point. 
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Fig. 8. — This curve relates perturbation size and inclination for the e Eri perturbation through the Pourbaix & Jorrisen (2000) relation 
(Equation 10). We use the curve as a 'prior' in a quasi-Bayesian sense. Our final values for the semimajor axis of the astrometric 
perturbation, 04 > and inclination, i are plotted with the formal errors. 
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Fig. 9. — Radial velocity measurements of e Eri from s our ces as indicated in the legend (and identified in Table 131 , The dashed line is 
the velocity predicted from the orbital parameters ( Table ITT1 derived in the combined astrometry and radial velocity solution (Section 6). 
The linear trend in velocity can be seen by comparing the minima in the radial velocity orbit near 1981 with that at 2001. 
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Fig. 10. — Perturbation orbit for e Eri. Elements are found in Table ITTI Residual vectors are plotted, connecting each normal point 
residual to its predicted position at each epoch of observation. Circle size is proportional to the number of observations forming the normal 
point. The square marks periastron passage, Tq=2007.3. The arrow indicates the direction of motion of e Eri. 



